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ABSTRACT 

A power-law density model, i.e., p{r) oc has been commonly employed in 
strong gravitational leasing studies, including the so-called time-delay technique used 
to infer the Hubble constant Hq. However, since the radial scale at which strong leas¬ 
ing features are formed corresponds to the transition from the dominance of baryonic 
matter to dark matter, there is no known reason why galaxies should follow a power 
law in density. The assumption of a power law artificially breaks the mass-sheet de¬ 
generacy, a well-known invariance transformation in gravitational leasing which affects 
the product of Hubble constant and time delay and can therefore cause a bias in the 
determination of Hq from the time-delay technique. In this paper, we use the Illustris 
hydrodynamical simulations to estimate the amplitude of this bias, and to understand 
how it is related to observational properties of galaxies. Investigating a large sample 
of Illustris galaxies that have velocity dispersion ctsie ^ 160kms“^ at redshifts below 
z = 1, we find that the bias on Hq introduced by the power-law assumption can reach 
20% — 50%, with a scatter of 10% — 30% (rms). However, we find that by selecting 
galaxies with an inferred power-law model slope close to isothermal, it is possible to 
reduce the bias on Hq to ^5%, and the scatter to <10%. This could potentially 
be used to form less biased statistical samples for Hq measurements in the upcoming 
large survey era. 

Key words: gravitational leasing: strong - galaxies: haloes - galaxies: structure - 
cosmology: theory - dark matter. 


1 INTRODUCTION 

Strong gravitational lensing is a major tool of mod¬ 
ern ext ragalactic astrophy sics. Ever since its first dis¬ 
covery JWalsh et al.l Il979l l. it has been used as a nat¬ 
ural telescope to magnify the distant Universe, a scale 
to weigh galaxies, and a ladder to measure the Hub- 
ble constant Hq - the expansion rate of the Universe 
llRefsdall Il964 l. With the current and upcoming space 
missions like GAIA and Euclid, and ground-based facili¬ 
ties such as the Large Synoptic Survey Telescope (LSST) 
and the Square Kilometre Array, strong lensing stud- 
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ies will experience an unprecedented opportunity to ex¬ 
ploit seve ral thousands of lensing galaxies to be discov¬ 
ered (e.g., Coe fc Moustakaj 20091: Qguri fc Marshall 20101: 


Finet et al.ll2012l : lLSST Dark Energy Science CoiTaborationl 

2012ll . while only hundreds are currently known. This tran¬ 
sition into the big data era demands a good understanding 
of and control over systematic errors that are introduced 
by lens modelling techniques, which will affect scrutinized 
quantities such as the lens density slopes and their evolu¬ 
tion with time, and derived cosmological parameters such 
as Ho- 


The observational properties of galaxy-scale strong lens 
systems (i.e. positions and flux ratios of unresolved images, 
brightness distribution of lensed extended components, and 
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time-delays) depend mainly on the mass distribution in¬ 
side and near the Einstein radius of the lens, which cor¬ 
responds in general to a few half-light radii. Within this 
radius, both dark and baryonic matter are believed to co¬ 
exist in roughly similar amounts. However, because of the 
interplay between cooling and heating mechanisms, the un¬ 
known efficiency of star formation, and a wide variety of 
feedback processes, a large scatter in the galaxy total mass 
profiles is naively expected. Observationally, results from 
galaxy dynamics combined with strong lensing measure¬ 
ments, but also X-ray emission of elliptical galaxies, sug¬ 
gest a limited diversity of profiles and indicate that the 
inner prohle is almost isothermal (e.g., iRusin e t al.l l2003l: 


HumDhre^_et_^ 2006|;^ Gnvazzi et al.ll200'J Koopmans et al.l 


20091: lAueer et al.ll2010ir This surprising result, also known 

as the “bulge-halo conspiracy”, sets impo rtant constraints 
on the formation history of galaxies (e.g., Koo pmans et al.l 


l2006l : ISonnenfeld et al]l201,9l : lDutton fc Treull2014h . In addi- 

tion, it has an important impact on numerous phenomeno¬ 
logical studies of galaxies, as it motivates simple power-law 
models i.e., p(r) oc r ' of the volume mass density pro¬ 
file, with a corresponding convergence profile k{9) oc 9^ ' 
[i.e. the surface density distribution S(9) normalized by the 
lensing critical density 

= ( 4 ^) (rZk) ’ 


where Dd, and Dds are the angular diameter distances 
to the lens, to the source, and from the lens to the source, 
respectively]. This model is commonly used in studies of 
galaxy-scale strong lens systems, and often isothermality, 
i.e., 7 ' = 2 , is assumed. 

Despite observational support, the power-law density 
profile is no more than a convenient “working model” 
for strong lensing studies. This family of density profiles 
as found by htting strong lensing measurements does not 
uniquely/necessarily represent the mass distribution of the 
inner regions of galaxies due to th e so-called mass-sheet de¬ 
generacy fMSD: iFalco et al.]ll985l ~l. Any given convergence 
field k{6) can be transformed into a new convergence prohle 
k\(9) via 


«:a( 0) = AK(0)-f (1-A). (1) 


Together with an isotropic rescaling of the source plane 
P A/3, the transformation (1) maintains the invariance of 
all observable quantities, such as image positions and flux ra¬ 
tios, except the time delays. This transformation is referred 
to as the “mass-sheet transformation” (MST). The only ob¬ 
servable quantity modihed under the MST is the product 
of the Hubble constant Hq and the time delay At, which is 
transformed such that HoAt —>■ XHoAt. Because of the in¬ 
variance of observables, information from gravitational lens 
properties cannot distinguish between any of the members 
of the family of mass proHles ka(0). However, if a particu¬ 
lar parametrized mass model is chosen, such as a (local or 
global) power-law prohle, the MST is artihcially broken by 
yielding that member of the family kx which is closest to the 
chosen mass prohle parametrization. This artihcial breaking 
of the MSD can lead to a systematic bias on the determi¬ 
nation of Ho, for example, if the logarithmic density prohle 
near the Einstein radius is systematically curved upwards or 
downwards relative to a power-law prohle, then the values of 


the transformat ion parameter A will be s ystematically above 
or below unity llSchneider fc Slusjl2013ll . 

The range of values of A expected from a population 
of len sing galaxies is currently unknown. ISchneider fc Slusd 
1120131) have presented an example of a realistic galaxy 
model, composed as the sum of a baryonic and a dark 
matter prohle, which could be transformed into an ap¬ 
proximate power-law density prohle through Eq. © with 
A ~ 1.2, leading to a syst ematic error of about 20% on Hq. 
Van de Ven et al. (l2009l ) have shown that a large variety 
of observationally motivated composite models of galaxies 
systematically produce prohles that are slightly more con¬ 
vex/concave than isothermal in regions probed by lensing. 
Here, we study the expected distribution of surface den¬ 
sity prohles of galaxies drawn from a cos mological hydrody¬ 
nami c simulati on - the Illustris Project (|Vogelsberger_gt_alJ 
2014I. see als o Vogelsberger et ^ l2014l : Genel et al.l 2014 ; 


Nelson et al.ll2015l ). By selecting galaxies that most resem¬ 


ble the observed lensing galaxies in term of mass/velocity 
dispersion, we test the validity of the power-law assump¬ 
tion in the central regions of galaxies, and evaluate the ex¬ 
pected distribution of the multiplicative bias A on Hq for this 
sample. We emphasize that in this analysis we assume that 
no other external constraint than the lensing information is 
available to break the MSD. Whereas stellar dynamical in¬ 
formation from spectroscopy can yield additional constraints 
on the lensing mass distribution and thus helps to limit the 
allowed range of A, the accuracy of this method is limited, 
e.g., by spectroscopic resolution and the unknown distribu¬ 
tion of stellar orbits. 

The outline of the paper is as follows. In Sect. 2, we 
describe the methodology to measure the slope and curva¬ 
ture of the surface density profiles of the simulated galaxies, 
and how we derive the A parameter which transforms the in¬ 
trinsic distribution into an approximate power law. In Sect. 
3, we give a brief description of the lensing galaxy samples 
that we have extracted from the Illustris simulation; general 
properties of the simulated galaxy sample are presented. As 
the (surface) density distributions are of particular interest, 
we dedicate Sect. 4 to a detailed discussion about the diver¬ 
sity of the central density prohles in regions probed by strong 
lensing. In Sect. 5, we present distributions of A that result 
from the use of power-law density prohles to model various 
samples of lensing galaxies. Such distributions also indicate 
by how much the derived value of Ho would be systemati¬ 
cally biased from the true value by employing the power-law 
assumption. A hnal discussion and conclusions are given in 
Sect. 6 . 


2 METHODOLOGY 

In this section we explain the quantities that we use to char¬ 
acterize the density prohles of the simulated galaxies drawn 
from Illustris. Since the MST only affects the radial density 
prohle, but leaves the shape of the isodensity contours invari¬ 
ant, we are only concerned with the azimuthally-averaged, 
i.e., radial density prohle k{9) of a lensing galaxy. 

In order to quantify how well power-law models can de¬ 
scribe the surface density prohle of a lensing galaxy, we have 
hrst calculated a mean local slope s and the corresponding 
curvature ^ following the treatment of Schneider & Sluse 
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(2013). Assuming that the strong lensing features (multiple 
compact images, images of extended source components) are 
located between an inner and an outer radius of 6 i = xi 0 b 
and 02 = X 2 dE, respectively, where 0e is the Einstein radius 
of the lens, a mean logarithmic slope s for the convergence 
distribution k{ 6 ) between 6 i and 62 can then be calculated 
asO 


_ \n{K2/Kl) 
\n{ 6 i/ 62 )' 


( 2 ) 


where k\ = k( 0 i) and K 2 = ^{ 82 ) are the convergence at 
the inner and outer image positions, s = 1 corresponds to 
the projected slope expected for an isothermal distribution; 
s > 1 (s < 1 ) indicates the proHle steeper (shallower) than 
isothermal. Unless mentioned otherwise, we set xi = 0.5 and 
X 2 = 1.5 in this paper, which are characteristic values for 
observed lensed systems. 

We define the curvature parameter ^ as 





y/l<ilK2 


(3) 


which indicates the closeness of the profile k{9) between 9i 
and 62 to a power-law distribution: an exact power-law pro¬ 
Hle will have 5=1; while 5 > 1 (5 < 1) means that the pro¬ 
jected density profile is concave-upward (convex-downward), 
which lies above (below) the power-law interpolation be¬ 
tween 01 and 02 - 

The MST (1) transforms the original k{0) into a new 
distribution kx{0), with slope 


_ \n{Kx{02)/K\{0l)) 
“ ln( 0 i/ 02 ) 

and curvature 


5a 


li\{\/0102) 

\/Kx{0l)Kxi02) 


leaving image astrometry and flux ratios unchanged. Here 
and below, we specifically denote A as the value of the trans¬ 
formation parameter that transforms k{ 0 ) to a new distri¬ 
bution Kx{ 0 ), which satisfies the criterion that 5 a = 1 ; this 
condition yields 


^ _ K2 + Kl - 2^^K2Ki 

K2 -I- Ki — 25-^K2 «:i + (5^ — l)«:2«;i 

Let r = 0 \I 02 < 1, so that K 2 /K 1 = r® (from Eq. [^. Then 
by dividing the numerator and denominator of Eq. 11 by 
Aci, we can re-write A as a function of r, s, 5 and K 2 : 


l + r‘’- 25 r'^/" 

1 -I- J.S _ 2 5 fs/2 _j_ ^^2 _ 


(5) 


Note that in the case where A = 1,5 = 5a = 1- 
Here and in the following, we use 5 a = 1 as our power- 
law criterion: that an MST with A given by Eq. 
transforms k( 0 ) into an approximate power law kx{0) be¬ 
tween 01 and 02 is in such a sense that the three points 
(In 61 , ln«:A(^i)), (In 02 , lnKA( 02 )) and their mid-point 
(In x/ 0102 , In ka(a/ 0102 )) he on a straight line in the log¬ 
arithmic k( 0 ) plot. 

A meaningful MST requires that (i) the transformed 


^ The quantity that we have denoted here as s is actually the 
same as s of Schneider &: Sluse (2013). 


density profile kx{0 ) remains monotonically decreasing, 
which yields A > 0; and (ii) kx{0) is non-negative over 
the range considered, i.e., iix{ 02 ) > 0, which yields A < 
1/(1 — K 2 ). Note that typically 85% of the galaxies in our 
lens samples satisfy condition (i) and all satisfy condition 
(ii). See Table 1 for the exact fractions. 

The slope sx of the transformed profile, which by con¬ 
struction is an approximate power law between 0i and 02 
(i.e., 5 a = 1 ), can explicitly be related to the properties of 
the original profile. Indeed, by definition it obeys the rela¬ 
tion 

sx ^a(02) \K 2 T (1 — A) ^2 -t- (1 — '^)/'^ /^\ 

-:^- A.. + (l-A) = ., + (l-A)/A ’ 

From Eq. @ we obtain 


1-A 


(5^ — 1 )kiK 2 


(7) 


A Kl -I- K 2 — 2^-^KiK2 
Combining Eq. 0 with Eq. then yields 

_ K2 — 2^fi2^yK,lK2 -I- 5^«:iK2 
— 2^Ki^/kiK^ -h ^^KiK2 

By dividing the numerator and denominator of Eq. ® by 
K 1 K 2 we obtain 

s/2 ^ 


(8) 


r‘ - 2^r‘'^ -h 5 ^ 

1 ^ — 25 ^—-f 5 ^ 


5-r^ 

1 — 5 r'’/^ 


(9) 


1 + 


^r-s/2 _ 


Hence, the power-law slope after MST is 


sa = s -h In 





/ Inr . 


( 10 ) 


Thus, Sx is a function solely of the original slope s and the 
curvature parameter 5 - H 5 > li then (since > 1 and 
Inr < 0 ) sa < s, i.e., in this case the transformed profile is 
flatter than the original. This case also corresponds to A < 1, 
i.e., the mass sheet added has positive convergence, leading 
to a flatter profile after MST. 

It is worth noting that a different radial range [0i, 02 ] 
would yield different kx and A. This is simply because they 
depend on the slope s and curvature 5 which are defined in 
terms of the angular range where strong lensing measure¬ 
ments are available. As stated before, we focus on the range 
01 = 0.5 0E and 02 = 1.5 0e. Results for a different angular 
range (with 0 i = 0.8 0e and 02 = 1.2 0e) are also presented 
in Appendix]^ As can be seen, although measurements in 
these two cases are not exactly the same on a one-to-one 
basis, the statistical distributions of A are independent of 
the choice of the 0 i {i = 1 , 2 ). 

Equations (1) to (10) relate to the local convergence 
profile k{9). They were derived and discussed in the frame¬ 
work of the assumption that the convergence profile follows 
approximately a power law in the range 0i ^ 0 ^ 02 . Al¬ 
ternatively, we can also define analogous quantities under 
the assumption that the radial dependence of the deflec¬ 
tion angle behaves like an approximate power law, which is 
equivalent to a power law in the mean convergence within 
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e, i.e., 

= e'Kie'). 

Hence, we can consider an MST that results in a transformed 
cumulative distribution Rx to satisfy our power-law criterion 
so that the profile of R\ is a good approximation of a power 
law between 6 i and 62 - The MST for k is the same as that 
for K, i.e., 


where baryons yield a significant contribution to the to¬ 
tal mass profile. On the other hand, the Illustris simulation 
provides a large and realistic sample of galaxies as needed 
to quantify in a statistical way the systematics associated 
with the assumption of power law density profiles in strong- 
lensing studies. 

3.1 Main samples of lensing galaxies 


ka = Ak -I- (1 — A). (11) 

Note that the radii where k = kx — 1 and where R = Rx ~ ^ 
remain unchanged; the latter corresponds to the Einstein ra¬ 
dius by definition. Hereafter we denote by barred symbols 
all quantities that are associated with the cumulative distri¬ 
bution R and its MST. The whole formalism from Eq. (2) to 
(10) also holds for R after replacing k with R, s and ^ with 
s and f, sx and ^x with sx and ^x, and A with A. 

Which of the two MSTs - the one that yields an approx¬ 
imate local power law of the convergence, i.e., = 1 , or that 

leading to an approximate power law in the deflection, i.e., 
^A = 1 “ is the more relevant one depends on which assump¬ 
tions are made in lens modelling. In many applications, it is 
assumed that the density profile follows a global power law 
(at least up to the radius where strong lensing features are 
observed), in which case the requirement ^a = 1 applies. Al¬ 
ternatively, lens models can be considered in which the local 
profile follows a power law, but an additional mass compo¬ 
nent must be assumed on angular scales smaller than 9i, i.e., 
the presence of a supermassive black hole, or a finite core 
radius of the mass distribution, both of which cause devia¬ 
tions from a global power law. Note that the assumption of 
a global power law implies that the density profile between 
9i and 02 is a local power law (with the same slope). Thus, 
the global power-law assumption is a stronger one than just 
enforcing a local power law. However, since the density pro¬ 
files of real galaxies do not follow power laws, this conclusion 
does not hold for real lenses. We note that in general, A 7 ^ A, 
and sx^sx- 


3 THE ILLUSTRIS LENS SAMPLES AND 

THEIR GENERAL PROPERTIES 

The Illustris Project is a series of cos mological hydrody- 
nami cal simula tions of galaxy formation ([Vogelsbergeret^^jj 
20 14 see al so IVqgelsberger et al. I |2014 iGenel et al.l 12014 
Nelson et al. ll 2 Ql^. The highest resolution run covers a vol¬ 
ume of (106.5 Mpc)® and has a dark matter mass resolution 
of 6.26 X 10® Mq and an initial baryonic mass resolution of 
1.26 X 10® M 0 , resolving gravitational dynamics down to a 
physical scale of e = 710 pc. Taking into accoun t various 
baryonic processes (IVoeelsberger et al.ll20l4 120l4l . such as 
gas cooling, stellar evolution and feedback, chemical enrich¬ 
ment, supermassive black hole growth and feedback from 
active galactic nuclei, the Illustris simulation resolves 40 000 
galaxies with a variety of morphologies and reproduces many 
fundamental properties of observed galaxies (e.g., galaxy lu¬ 
minosity functions and Tully-Fisher relations). 

The Illustris Project is an ideal “laboratory” for our 
study. On the one hand, this is because galaxy-scale strong 
lensing probes the projected central regions of galaxies 


We take the lensing galaxies at various redshifts from 
the mock strong-le ns catalogue of the Illustris simulation 
l|Nelson et al.ll2()la . see also Xu et al., in preparation). Our 
selection criteria are as follows: for a given source redshift Zs, 
the angular (and physical) Einstein radius 0 b (and Re) of 
each galaxy is first determined as the radius within which the 
mean convergence R = 1. Referring to the singular isother¬ 
mal sphere model, we characterize each lens galaxy by its 
“velocity dispersion” usie, which is linked to the Einstein 
radius via 0 b = 47 r(crsiE/c)^(Dds/Ds). We select galaxies 
that have ctsib ^ 160kms“^. This lower limit on the veloc¬ 
ity dispersion is motivated by the observed lens sample from 
the SLAGS survey, where stellar velocity dispe rsions of lens 
galax ies range from 160 km s”*^ to 400kms“^ (iBolton et al.l 

l2008ll . 

The Illustris simulation does not sample the very mas¬ 
sive end of the mass spectrum well, i.e., galaxy clusters. 
There are a couple of massive systems that have total masses 
around We have adopted a further selection 

criterion, i.e., excluding all satellite/companion galaxies in 
the group and galaxy environment and only selecting the 
“central” galaxies to form our sample. The projection ef¬ 
fects from their companions (as well as galaxies along the 
line of sight) have been excluded from the calculation of the 
density profiles of these “central” galaxies. 

Note that the surface brightness distributions of the Il¬ 
lustris lensing galaxies can be well fitted by Sersic profiles 
llSersiclll963l l. In particular for the elliptical galaxiefl whose 
stellar masses and velocity dispersions are in the same ranges 
as the SLAGS samples, their logarithmic density slopes are 
also consi stent with the one constrained for the SLAG S sam¬ 
ples fe.g.. lAueer et al.]|2010l : ISonnenfeld et al.|[2013l ~l. A de¬ 
tailed comparison of Illustris galaxies to the SLAGS sam¬ 
ple will be presented in a forthcoming paper (Xu et al., in 
preparation). In this paper, as we are interested in the sur¬ 
face density distribution, which is a projected quantity, and 
in order to increase the sample size, we treat the three inde¬ 
pendent projections of each selected galaxy as independent 
lenses. This yields more than 1000 surface density profiles 
for the analysis at each of the lens redshifts studied. 

In the following we refer to the region where strong 
lensing measurements are available (i.e., the range [ 0 i, 02 ]) as 
the strong lensing region. Because 0 b is an angular quantity, 
the strong lensing region depends on the lens redshift Zd 
and the source redshift 2 s. In order to separate the effect 
of galaxy evolution from a change of Re (such that lensing 
probes different regions of a same galaxy), we have formed 


^ The surface brightness distribution of each galaxy in the mock 
catalogue has been fitted using both the de Vaucouleurs and the 
exponential disc profiles. A galaxy is classified as an early (late) 
type if the former (latter) provides a better fit. 



























Lens profiles, power-law models and bias on Hq 5 


Table 1. A summary of sample properties (galaxies in their three independent projections with ctsie ^ 160kms“^ are selected). The 
total number of galaxy projections that meet the selection criterion are given in row (3); row (4), (5) and (6) give the percentages of 
galaxies whose profiles can be transformed to power laws via MST (see Sect. 2); in row (7) is the minimal physical Einstein radius, 

which corresponds to the lower limit of o-giE = 160kms“^; in row (8), (9) and (10), the mean and medi an Einstein radius and its 
standard deviation ctr^ are presented, respectively. 


Sample sets 


-2s 

= 1.5 


2d = 0.6 

Redshifts 

2;J = 0.2 

O 

II 

2d = 0.6 

2d = 0.8 

2s = 1.0 

Zs — 1.5 

2s = 3.0 

Total number of projections 

1044 

1334 

1433 

1363 

1092 

1433 

1673 

Meaningful MST for k 

90% 

93% 

90% 

88% 

88% 

90% 

93% 

Meaningful MST for R 

77% 

82% 

78% 

75% 

73% 

78% 

81% 

Meaningful MST for both k and R 

76% 

81% 

77% 

73% 

72% 

77% 

80% 

(kpc) 

2.00 

2.60 

2.50 

2.06 

1.68 

2.50 

3.29 

Mean Re (hpc) 

4.16 

5.32 

5.14 

4.41 

3.71 

5.14 

6.47 

Median Re (kpc) 

3.62 

4.69 

4.49 

3.91 

3.29 

4.49 

5.60 

Standard deviation (kpc) 

2.04 

2.53 

2.46 

2.04 

1.75 

2.46 

3.15 


two sets of samples according to different combinations of 
Zs and Zd- 

In the first set, we have selected galaxies at Zd = 
[0.2, 0.4, 0.6, 0.8], assuming a fixed source redshift at 
Zs = 1.5 (Tabled] left-hand columns). Note that assuming 
the innermost boundary of the strong lensing region to be 
0.5i?E, the Zd — Zs combinations used here guarantee that 
even the smallest galaxies (with Einstein radius of 
which corresponds to the lower limit of ctsie = 160kms“^) 
are fully resolved at radii larger than the innermost bound¬ 
ary, i.e., 0.5i?E‘“>e (where e = 710 pc is the simulation 
softening length). As can be seen from Table dl the mean 
Einstein radius of the galaxies, (Re), is approximately the 
same for the four Zd — Za combinations. We therefore probe 
with this sample typically the inner 5kpc of galaxies, such 
that any difference in the final statistical results is due to 
galaxy evolution. 

In the second sample set, we fix the lens redshift at 
Zd = 0.6 and assume source redshifts at Zs = [1.0, 1.5, 3.0] 
(Table dl right-hand columns). Again the Zd — Za combi¬ 
nations applied in this sample set also guarantee that the 
smallest galaxies are fully resolved at radii larger than the 
innermost boundary of their strong lensing regions. For fixed 
Zd, the higher Re correspond to the larger Za'. the mean (Re) 
for the Za = 3.0 sample is nearly twice as large as that for the 
Za = 1.0 sample. This allows us to study the dependence of 
our results on the radii probed by strong lensing. Note that 
different numbers of galaxies are obtained in spite of the 
same ctsie cut, because the actual galaxy mass distributions 
are not truly isothermal, which thus introduces scatter in the 
crgiE-mass relation. The scatter in mass that corresponds to 
CTSIE = 160kms“^ for the 2s = 1-0 sample results in a mini¬ 
mum halo mass of ~ 9x Mq-, while the scatter in the 

2 s = 3.0 sample leads to a slightly lower halo-mass limit of 
~ 7 X 10^^ The latter sample therefore has a larger 

number of galaxies than the former. 


3.2 Galaxy profiles and general properties 

In the following, we present the surface density profiles of 
the selected galaxies and some of their general properties. 



Figure 1. The s — ^ distribution of the selected lensing galaxies 
at Zd = 0.2, assuming Zs = 1.5. The sample is divided into four 
subgroups according to their s and ^ values: subgroup I (red) is 
defined as galaxies with s > 1.2, ^ > 1.02; subgroup II (blue) 
are those with s > 1.2, ^ < 0.98; subgroup III (orange) repre¬ 
sents galaxies with s < 0.95, ^ < 0.98; and subgroup IV (green) 
are those with s < 0.95, ^ > 1.02. Contours indicate where the 
transformed slopes sx = [0.3, 0.5, 0.7, 0.9, 1.1, 1.3, 1.5, 1.7]. 

For each galaxy, we have azimuthally averagecQ the surface 
mass density distribution to obtain the radial profile of the 
convergence in the projected central region. To derive the 
profile parameters, we fit a 10-order polynomial to the log- 


^ We have also calculated the radial profile by averaging the con¬ 
vergence distribution in elliptical annuli, the shape of which is 
determined by the second moment of total projected mass within 
20 e- We have verified that both the curvature parameter ^ and 
the MST parameters A under the elliptical average are consistent 
with the ones derived under azimuthal averaging; the final result 
remains. 
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arithmic radial profile within a radius range from 2e to five 
times the half-stellar-mass radius of the galaxy, which well 
covers the region constrained by strong lensing. Using the 
polynomial htted radial profiles, we measure the slopes and 
curvatures of both the local convergence k and the cumula¬ 
tive distribution R, which are then used to derive the MST 
parameters A and A and the transformed slopes s\ and s\, 
according to the formalisms derived in Sect. 2. 

Figure [1] shows the distribution of s versus ^ for the 
lens sample with za = 0.2 and = 1.5. The overlaid con¬ 
tours indicate the transformed slope s\ as a function of (s, 
^), which increases from the bottom right quadrant to the 
top-left quadrant in the s — ^ plane. The mean slope s mea¬ 
sured between 0.5 0b and 1.5 0e typically ranges from 0.5 to 
1 .6, while for the curvature we have |C ~ 1| 15% for all 

our samples. Of course, the range of ^ strongly depends on 
A0 = 02 — Oi: as A6 —>■ 0, ^ ^ 1. Note that galaxies that 
fail the criteria “A > 0” (and thus sa > 0) for a meaning¬ 
ful MST lie at the lower-right corner of the diagram, to the 
right of the green points. These galaxies are not fundamen¬ 
tally different from the “green” galaxies in terms of their 
physical properties except that they have even shallower 
profiles with larger curvatures. Mathematically in order to 
transform their profiles into power laws, the additional mass 
sheet (1-A) has to be so large that the transformed power- 
law densities become lower in the centres and larger at the 
outskirts, and thus sa becomes negative. 

In order to understand different profile behaviours, we 
have further divided each lens sample into four representa¬ 
tive subgroups according to their s and ^ values. As shown 
by the coloured symbols in Fig. [T] subgroup I as indicated 
by red is defined for galaxies with s > 1.2 and ^ > 1.02, 
i.e., galaxies with mean slopes steeper than isothermal and 
having concave-upward profiles; subgroup II as shown in 
blue consists of galaxies with s > 1.2 and ^ < 0.98, i.e., 
galaxies that also have mean slopes steeper than isothermal 
but convex-downward profiles; subgroup III as marked in or¬ 
ange represents galaxies with s < 0.95, ^ < 0.98, i.e., those 
with mean slopes flatter than isothermal and with convex- 
downward profiles; and subgroup IV plotted in green is com¬ 
posed of galaxies with s < 0.95, ^ > 1.02, i.e., those with 
mean slopes flatter than isothermal but having concave- 
upward profiles. 

Since ^ = 1 indicates approximate power-law profiles, 
it is used to separate subgroups. The division using s = 1 is 
somewhat arbitrary; this is purely motivated by the obser- 


vation that the distribution of the total density s 
galaxies seems to peak at the isothermal slope ( 

lopes of lens 

Rusin et al.l 

20031: iRusin & Kochanekll2005l:lKooDmans et al 

l2006l.l2009l: 

Gavazzi et al.l 20071: Auger et al.l 2010l). In other words, the 


four subgroups are extreme cases where the density profiles 
are far from isothermal and perfect power law. 


Apart from the profile parameters, for each galaxy in 
projection we have also calculated a variety of properties: 
(1) the total stellar mass M* of the galaxy and the total 
mass Mtot of the galaxy plus its host dark matter halo; (2) 
the effective radius 0eff, within which half of the galaxy lumi¬ 
nosity is enclosed; (3) the cumulative dark matter fraction 
/dm(^ 0 e) projected within 0b and the local dark matter 
fraction /dm(0E) projected at9^\ (4) the angular radius 0/5o 
where the projected density distribution of dark matter in- 
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Figure 2. The distributions of Mtot versus ctsie (top) and /dm{^: 
0e) versus Mtot (bottom). The galaxy sample as well as the colour 
coding is the same as used in Fig. [T] 


tercepts that of the baryonic matter, i.e., the radius where 
the local dark matter fraction /dm(0/5o) = 0.5. 

Fig.[2]shows the distributions of Mtot versus ctsie (top) 
and /dm(^ 0 b) versus Mtot (bottom) for the same galaxy 
sample as in Fig.[T] Different projections of the same galaxy 
result in multiple values of ctsie and /dm(^ 0b) at the same 
mass, reflecting the triaxiality of the lens mass distribution. 
The cumulative dark matter fraction /dm(^ 0b) typically 
ranges from 30% to 70% with more massive galaxies hav¬ 
ing higher /dm(^ 0 e), in good agreement with strong lens¬ 
ing observati ons (e.g ., |Treu_J^J4oopmaJij 2004 ; Koopmans 
et al. 2006: Jimg&jCochanek 20071: Tortora et al. l2009l: 
Gjdllo et al.l 2003 : ICardone et al. 20091 : Nanolitano et all 

201ol). 

It can be seen clearly that the four subgroups have 
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different galaxy properties, e.g., the shallower- and convex- 
profiled galaxies (subgroup III) are among the most massive 
ones and have the largest velocity dispersions and the high¬ 
est dark matter fractions. In the next section, we aim at 
understanding the direct reason for such a diversity in the 
surface density profiles measured in strong leasing regions. 


4 THE CENTRAL SURFACE DENSITY 

PROFILE: SLOPE AND CURVATURE 

Galaxy-scale strong leasing probes the central few kpc re¬ 
gion of a leasing galaxy. The central (surface) density profile 
is of particular interest, because on such scales, the enclosed 
mass of dark and baryonic matter are quite similar. Various 
processes associated with star formation and feedback to¬ 
gether shape the density distributions of both components 
and their sum. 

We first illustrate in Fig. [3] the projected density pro¬ 
files associated with the four galaxy subgroups fSect. |3^ . 
For each subgroup we show ten typical density profiles and 
their dark matter fraction distributions. The projected radii 
are normalized to 0 b (see Fig. IB2I for profiles as a function 
of 0eff). In the top panels, the black curves are the total 
surface density distributions k{9), scaled up by a factor of 

5 for clarity; the red and blue curves represent profiles of 
the projected dark matter K( 0 )/dm( 0 ) and projected bary¬ 
onic mass ft(0)[l — /dm(0)], respectively; the dashed green 
line indicates the logarithmic slope of s = 1 from 0.5 0 e to 
1.5 0E. In the bottom panels, the cumulative and local dark 
matter fraction /dm(^ 0) and /dm(0) are given in black and 
green, respectively. 

The first conclusion is that, as expected, the baryonic 
distribution is much more concentrated than the dark mat¬ 
ter distribution. In regions probed by strong lensing, the 
two profiles intercept; the slope of the former is in general 
steeper than that of the latter. As a result, baryons and dark 
matter dominate different parts of the total (surface) den¬ 
sity profile: the former in the inner regions and the latter 
at larger radii; the dark matter fraction therefore increases 
with radius. 


4.1 Individual profiles of the two ingredients 

As the total density profile is the sum of two components, we 
present hereafter slopes s and curvatures ^ of the baryonic 
(i.e. Sb and ^b) and dark matter profile (i.e. Sdm and ^dm) 
separately. 

Fig. [4] and [5] show the dependencies of the slopes and 
curvatures on 0/5o/0e and /dm (0 b), which are two key quan¬ 
tities that are closely related to different profiles of sub¬ 
groups. In particular, 0/5o/0e describes where, with respect 
to 0 B, the fraction of the projected dark matter catches up 
and the distribution intercepts that of baryons. Different val¬ 
ues of 0/50/0B reflect different parts of the profile that strong 
lensing probes: a smaller 0/5o/0b means that the strong lens¬ 
ing region (~ 0 b) is at larger radii than the baryon-dark 
matter interception radius 0 / 5 o, and thus the dark matter 
fraction /dm (0 b) is larger. As indicated by the colour cod¬ 
ing, in the four highlighted subgroups, strong lensing probes 
different parts of the density profile, from much further-out 


in subgroup III (orange) to most closer-in in subgroup IV 
(green). 

As can be seen from Fig. 21 the central density pro¬ 
file of the dark matter distribution is always shallower than 
isothermal, i.e., Sdm < 1 ; while that of the baryonic mat¬ 
ter in most cases is considerably steeper, i.e., Sb > 1- From 
subgroup IV (green) to I (red), to II (blue) and to III (or¬ 
ange), as strong lensing probes increasingly larger radii of 
the galaxy, a significant profile steepening followed by a slope 
flattening develops with increasing /dm (0b) (or decreasing 
0/50 /0e); this is seen for both dark matter and baryonic 
components in similar fashion. 

In order to study the correlation between dark matter 
and baryonic slopes in the strong lensing region, Fig.[ 6 ]shows 
the projected baryonic slope Sb as a function of dark matter 
slope Sdm for the galaxy sample with Zd = 0.2 and Zb = 1.5. 
In fact a strong correlation between the two slopes is found 
for all of our samples with different Zd — Zs combinations. 
This is a consequence of the strong interplay between dark 
matter and baryons in the central regions of galaxies. 

The dependencies of the curvature parameters, as pre¬ 
sented in Fig. [5] also exhibit similar distributions between 
dark matter and baryons. Once again, as the strong lens¬ 
ing region “moves out” (i.e., 0/5o/0e decreasing), a slight 
increase (with concave curvature) followed by a decreas¬ 
ing trend (to becoming convex) develops from subgroup IV 
(green) to I (red), to II (blue) and to III (orange). 

Figures[4]and[5]indicate the general trend of the (central 
and logarithmic) surface density profiles of both dark matter 
and baryonic distributions: with an increasing radius from 
the centre of a galaxy, a shallower and concave-upward pro¬ 
file in the inner region (probed by strong lensing of subgroup 
IV in green) is followed by a much steeper profile at inter¬ 
mediate radii (now probed by strong lensing of subgroup 
I in red and II in blue), and eventually is turning into a 
(slightly) shallower profile again with convex curvature at 
larger radii (now probed by strong lensing of subgroup III 
in green). Strong lensing therefore probes different regions 
of such matter density distributions. In a way, the distinc¬ 
tive profiles shown in Fig. Oof galaxies in the four different 
subgroups are actually “snapshots” taken at different parts 
of the density distributions. 

4.2 Total density profiles 

We next present the results of the total (surface) density pro¬ 
files. Fig. [7] shows the distributions of s versus 0/5o/0b (top 
left), s versus /dm (0b) (top right), ^ versus Ofso/ds (bottom 
left) and ^ versus /dm(0E) (bottom right). Compared with 
Fig. [4] and [5] the distributions of the total density profile pa¬ 
rameters have maintained the trends that are seen for both 
the dark matter and baryonic components. 

Strong gravitational lensing probes different parts of the 
total density profile. This is the direct reason why we see 
a collection of various profile behaviours. Among the four 
subgroups, the most massive subgroup - the one with shal¬ 
lower and convex density profiles (III, in orange) - have the 
smallest 0/5o/0e. In particular their 0/5o.$ 0.50 b indicates 
that their dark matter-baryon interception radii are located 
within the inner boundaries of the strong lensing regions. 
As a result, the total density profile measured within the 
strong lensing regions is strongly dominated by dark matter. 
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Figure 3. From each subgroup of the lens sample used in Fig. ^ ten typical surface density profiles and projected dark matter fraction 
distributions are shown. In the top panels the black curves show the total surface density distribution scaled up by a factor of 5 for 

clarity; the red and blue curves represent profiles of the projected dark matter «^(^)/dm(^) projected baryonic mass k,{ 0)[1 — fdmW]^ 
respectively; the dashed green line indicates the logarithmic slope of s = 1 from O.50 e to 1.50 e- the bottom panels, the distributions 
of the cumulative and local dark matter fraction /dm{^ /dm(^) ^re given in black and green, respectively. 
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Figure 4. Top left: dark matter slope Sdm versus ^/so/^E- Top 
right: Sdm versus /dm(^E)- Bottom left: baryonic slope Sb versus 
^/ 5 o/^E- Bottom right: Sb versus /dm(^E)- The galaxy sample as 
well as the colour coding is the same as in Fig. ^ 


Figure 5. Top left: dark matter curvature ^dm versus 0/5o/^E- 
Top right: ^dm versus /dm(^E)- Bottom left: baryonic curvature 
iJb versus ^/sq/^e- Bottom right: ^b versus /dm(^E)- The galaxy 
sample as well as the colour coding is the same as in Fig. ^ 


whose distribution is considerably shallower than isothermal 
(see Fig. mi. In comparison, subgroup II galaxies (in blue) 
have 0/50 > 0.50b, i.e., the dark matter-baryon interception 


radii now exceed the inner radii of the strong lensing re¬ 
gion, so that baryons can dominate the total density dis¬ 
tribution at 01 (see Fig. 3), resulting in the total slope s 
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Figure 6. The projected baryonic slope st, versus dark matter 
slope Sdmi measured between d\ = 0.5 0b and 02 = 1.5 0b. The 
galaxy sample and the colour coding are the same as in Fig. [J 
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steeper than isothermal. For both subgroups, a large frac¬ 
tion of the measured density distribution is dominated by 
dark matter; however, baryons contribute to a marked den¬ 
sity upturn around the inner boundary of the strong lensing 
region, resulting in a convex curvature < 1 ) of the total 
profile. 

In comparison, subgroup I (in red) and IV (in green) 
have even lower /dm( 0 E) and larger 0 / 6 o/fe, which means 
that baryons heavily dominate the strong lensing region in 
these cases. The behaviour of the total density profile, there¬ 
fore, largely depends on how the baryonic matter is dis¬ 
tributed. The strong lensing regions of subgroup IV galax¬ 
ies are much closer to (the projected) galaxy centres than 
their subgroup I counterparts. This can be seen from larger 
dfso/Gs of the former than of the latter. As shown in Sect. 
4.1, both dark and baryonic matter distributions are much 
shallower in the inner regions than at larger radii (see Fig. 
0 . As a result, the total density prohles of subgroup IV are 
shallower than of subgroup I galaxies; while both have con¬ 
cave curvatures in the strong lensing regions (i.e., ^dm > 1 
and 5b > 1, see Fig. 0 . 


5 THE DISTRIBUTIONS OF A AND THE 

CONSEQUENCE FOR Ho MEASUREMENTS 

We have shown in Sect. 0 and 0 that galaxies can have a 
significant concavity/convexity compared to power law den¬ 
sity profiles. Because of the existence of the MST, the lens 
modelling of these galaxies using a power-law density profile 
will yield a multiplicative bias A on Hq. We study hereafter 
the distributions of A resulting from the MST for the vari¬ 
ous samples of Illustris galaxies. Following our notation, the 
MST applied to the density profile k (via Eq. 0 such that 
the curvature parameter becomes unity, i.e., 5 a = 1 , yields 


Figure 7. Top panels: s versus 0/5o/0B (left) and s versus 
/dm(^B) (right). Bottom panels: 5 versus 0/5 o/0b (left) and 5 
versus /din(®E) (right). The galaxy sample and the colour coding 
are the same as in Fig. [T] 


a bias A and a transformed profile k\ with a slope sa, while 
the MST on k (via Eq. Illll is characterized by A and sa. 

In Fig. [S]we present three sets of distributions, i.e., A 
and A as a function of ctsib (left), 0B/^efr (middle), and as a 
function of /dm(41 ^b) (right), for the same galaxy sample as 
used before (where Zd = 0.2 and Zs = 1.5). We see that the 
distributions of A and A span a wide range of values from 0 to 
1.5 and exhibit a large scatter without a strong dependence 
on any of the observables above. Figures displayed in Ap- 
pendix[C]show that similar distributions are also present for 
all other studied 2 d — 2 s combinations. It is however notice¬ 
able that for galaxies with 200kms“^ < ctsib < 300kms“^, 
the mean multiplicative bias (A) is generally close to unity 
with deviations not larger than ~ 20 % at lower ctsib, while 
for A the mean deviation can be as large as ~ 50%. In all 
cases, an rms of 10% — 30% is present. This emphasizes 
that the power-law assumption in lens modelling is highly 
questionable, as it will create non-negligible biases in the 
derived values of Ho- However, we note that the observed 
“time delay lenses” have generally large velocity dispersion 
(ctsib ^ 250 kms“^), for which A is close to unity on average, 
though with a 10 % — 20 % scatter. 

Interestingly, when correlating A with sx (and A with 
Sa), which is the “measured” mean slope using power-law 
lens models, we see a much clearer systematic trend with 
smaller scatters. This is shown in Fig. 0 where the distri¬ 
butions of A versus s\ (left) and A versus s\ (right) are 
presented. Such a tight correlation is also seen for all other 
samples with different 2 d — Zg combinations (see Fig. EU. 

A striking result from Fig. 0 is that the s\ (and s\) 
distribution of galaxies with A « 1 (and A ~ 1) peaks at 
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Figure 8. The distributions of A (top) and A (bottom) versus ctsie (left), versus (middle) and versus /dm(^ Sg) (right). The 

sample as well as the colour coding is the same as in Fig.[T] Similar distributions are also seen for all other studied — Zs combinations. 
The upper axis of the top-left panel indicates the physical Einstein radius. 


sx ~ 1 (and s\ ~ 1); and the A (and A) distribution of 
galaxies with sx ~ 1 (and sx ~ 1) peaks at A « 1 (and 
A « 1). This could be potentially used to select time-delay 
samples which minimize the impact of the MST on Hq under 
the power-law assumption. 

We have further investigated such a possibility. For each 
of our lens samples, we have defined three subsamples that 
satisfy (1) sx £ [1 —As, 1 + As]; ( 2 ) s\ G [1 —As, H-As] and 
(3) both sx G [1 — As, 1 -I- As] and sx G [1 — As, 1 -I- As], 
respectively. The slope span As is set to be 0.1, i.e., the 
“measured” mean slopes (under the power-law assumption) 
are very close to isothermal. Table [2] summarizes the statis¬ 
tical properties of A and A distributions for the subsamples. 
In all the cases using these subsamples, both the systematic 
deviation from unity ( < 5%) and the scatter ( < 10%) are 
much smaller than those for the general samples. 

The statistical results presented above indicate that the 
power-law assumption will introduce non-negligible biases 
in Ho measurements, even though adopting power-law lens 
models is a common practice. However, through a careful 
sample selection, e.g., selecting galaxies that have larger ve¬ 
locity dispersions (i.e., large angular separations of the mul¬ 
tiple images) or those which have the “measured” profiles 
close to isothermal under the power-law assumption, the 
use of simple power-law models may yield estimates of Hq 
which are biased only at the few percent level - compared to 


a much larger bias from the total sample studied here. This 
work also shed light on such a possibility in this direction. 


6 DISCUSSION AND CONCLUSIONS 

The central regions of massive galaxies (i.e. typically the in¬ 
ner 5-10 kpc) contain dark and baryonic matter in roughly 
similar amounts. Although both components follow differ¬ 
ent density distributions, it has been found that the total 
density profile in that region is well approximated by a sin¬ 
gle power law (e.g., Koopmans et al. 2006, Gavazzi et al, 
2007). This has motivated the use of the power-law model 
p oc r as a generic model for numerous galaxy-scale 
strong lensing studies. However, as emphasized by Schneider 
& Sluse (2013), a mass distribution n{9) can be transformed 
into an approximate power law by means of an MST, i.e. 
ii\{9) — \k{6 ) -b (1 — A), even if the true mass profile shows 
considerable curvature. Under a MST, none of the lensing 
observables are modified except for the product of the Hub¬ 
ble constant and the time delay Hq At, which is scaled by 
the same factor A as the one characterizing the MST. By 
assuming a power-law profile, one “artificially” breaks the 
MSD and derives a biased estimate of Hq. In addition, the 
use of a single quantity, i.e., the logarithmic slope 7 ', to 
characterize the density profile of galaxies may be mislead¬ 
ing, since the intrinsic density does not follow an exact power 
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Figure 9. The distributions of A versus s\ (left) and A versus s\ (right). The sample as well as the colour coding is the same as in Fig. 
[T] Similar distributions are also seen for all other studied — Zg combinations (given in the Appendix[Cll. The solid lines indicate where 
A = 1 and where A = 1; while the dashed lines indicate where s> = 1 and where s\ = 1. 



Table 2. A summary of the statistical properties of A and A for galaxy subsamples that satisfy sx S [0.9, 1.1] and/or G [0.9, 1.1]: 


Sample sets 


Za = 1.5 


Zd = 0.6 

Redshifts 

Zd = 0.2 

Zd = 0.4 

Zd = 0.8 

Za = 1.0 

Za = 1.5 

Za = 3.0 

Subsample I: sx G [0.9, 1.1] 

Number of galaxy projections 

142 

210 

174 

146 

223 

296 

Mean A 

1.00 

0.98 

0.97 

0.97 

0.97 

0.97 

Median A 

0.99 

0.97 

0.96 

0.95 

0.96 

0.97 

Standard deviation ax 

0.11 

0.09 

0.10 

0.12 

0.10 

0.08 

Subsample II: sx £ [0.9, 1.1] 

Number of galaxy projections 

110 

167 

148 

102 

184 

278 

Mean A 

0.98 

0.98 

0.97 

0.96 

0.97 

0.98 

Median A 

0.98 

0.98 

0.95 

0.96 

0.97 

0.97 

Standard deviation a^ 

0.08 

0.09 

0.11 

0.09 

0.10 

0.14 


Subsample III: sy G [0.9, 1.1] and sx € [0.9, 1.1] 


Number of galaxy projections 

19 

36 

21 

19 

37 

63 

Mean A 

1.04 

1.02 

1.04 

1.06 

1.01 

1.01 

Median A 

1.03 

1.02 

1.01 

1.04 

1.02 

1.01 

Standard deviation ax 

0.06 

0.06 

0.10 

0.08 

0.05 

0.07 

Mean A 

1.02 

1.02 

1.04 

1.06 

1.01 

1.02 

Median A 

1.01 

1.01 

1.02 

1.04 

1.01 

1.01 

Standard deviation ax 

0.05 

0.06 

0.12 

0.08 

0.07 

0.08 


law. Therefore, such a slope depends on the radius, and its 
comparison among galaxies depends on the range of galac- 
tocentric radii over which it is measured [note that Dutton 
& Treu 2014 suggests the use of a mass-weighted slope to 
alleviate the problem]. 


In this paper, we have used a statistical sample of mock 
leasing galaxies from the first high-resolution cosmological- 
scale hydrodynamic sim ulation - the Illustris Project 
llVogelsberger et al.l 20141. see al so Vogelsberger et aj]|2014l : 
iGenel et al.l 20141 : Nelson et al.l l2015l ) - to test the valid¬ 
ity of the power-law assumption for the (central) surface 
density distribution of strong gravitational lenses. The sim¬ 
ulation reproduces extremely well numerous observational 
properties of real galaxies. In particular, we find that the 


cumulative dark matter fraction /dm(^ ^b) of the simulated 
leasing galaxies ranges from 30% to 70%, in good agreement 
with observations (Treu & Koopmans 2004; Koopmans et 
al. 2006). The average isothermality of the profiles derived 
for the simulated galaxies is also in agreement with various 
observational studies (e.g., Rusin et al. 2003; Auger et al. 
2010; Sonnenfeld et al. 2013). A detailed analysis of the 3d 
density proHle in the inner regions of Illustris lenses will be 
presented in a forthcoming paper (Xu et ah, in preparation). 

From the mock Illustris lensing galaxies at each of the 
studied lens redshifts below 2 : = 1, we have taken a sample of 
about 400 plausible lens galaxies, which have central velocity 
dispersion a ^ 160kms“^ and are “observed” through their 
three independent projections, yielding a sample of more 
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than 1000 projected density profiles. The slopes and curva¬ 
tures of the projected density profiles in regions probed by 
strong lensing (typically between 0.5 and 1.5 angular Ein¬ 
stein radii 9e) have been measured. We have transformed 
these profiles into (local) power laws via an MST and derived 
the distribution of the multiplicative bias A on Hq implied 
by a power-law assumption of the density profile. 

Our main findings are as follows; 

• The total surface density proEle in the projected cen¬ 
tral regions of galaxies depends on radius, showing devia¬ 
tions from a power law following some generic trends. At 
smaller radii the proEles are shallow and concave upward. 
With increasing radius, the profiles gradually steepen, and 
eventually at larger radii they become marginally shallower 
and attain convex curvatures. The projected density dis¬ 
tributions of both dark and baryonic matter follow such a 
general trend, while the latter is more concentrated (and 
steeper) than the former. Good correlations exist between 
the dark matter and baryonic density slopes throughout the 
projected radii and redshifts studied in this work, indicat¬ 
ing the strong interplay between dark matter and baryons 
in central regions of galaxies. 

• In the region probed by strong lensing, the (projected) 
total density proEle deviates signihcantly (by up to 15%) 
from an exact power law, with a variety of slopes and cur¬ 
vatures. The main reason for the diverse measurements of 
the profile parameters is that strong lensing probes differ¬ 
ent parts of this density distribution; the measured slopes 
and curvatures depend on the relative contribution of dark 
and baryonic matter in regions where the measurements are 
made. In general dark matter dominates the strong lens¬ 
ing regions of massive galaxies, which therefore tend to 
have mean profiles (in the strong lensing region) shallower 
than isothermal. In contrast, lower-mass galaxies have non- 
negligible baryonic components within their strong lensing 
regions; the shapes of the mean total density proEles are 
thus largely determined by the baryonic distributions. 

• The curvature of the (logarithmic) density profile is con¬ 
siderable and produces a non-negligible bias in the measured 
Hq, if only lensing information is used. The statistical dis¬ 
tribution of A (and A) derived from the Illustris strong lens 
sample spans from 0 to 1.5 and exhibits large scatter without 
a strong dependence on observationally-constrained galaxy 
properties, such as ctsib, dE/deS and /dm(^ fe)- In partic¬ 
ular, the mean deviation of A (and A) from unity can be 
as large as 20% — 50% with a scatter of 10% — 30% (rms) 
for galaxies with ctsie = 200kms“^ — 300kms~^. We note, 
however, that galaxies with derived slopes (under the power- 
law assumption) close to isothermal exhibit much smaller 
systematic deviation of < 5% from unity with a scatter of 
< 10% (rms). This could potentially be used to form a reli¬ 
able sample of lensing galaxies for Hq measurements in the 
upcoming big data era. 

We would like to stress that by finding an MST which 
approximates k\ as closely as possible by a power law, we 
in essence approximate k locally by a power law plus a 
uniform mass sheet. However, this is purely a mathemat¬ 
ical description; this uniform mass sheet is not assigned 
any physical meaning. In particular, this mass sheet can¬ 
not (or only partly) be interpreted as a convergence due 
to foreground and background material along the line of 


sight. The latter contribution to the convergence can be es¬ 
timated/constrained using cosmological simulations and/or 
through det ailed studies of the len s environment or weak 
lensing fe.g.. lSuvu et al.ir2010l . 1201^ 1. Furthermore, even if a 
lens is perfectly fitted with k being a power law plus uniform 
mass sheet, then the same is true for all k\ through MST. 
Hence, A remains undetermined. 

It is also worth noting that one may try out different 
“realistic” model profiles for the mass distribution of l enses , 
such as the comp o site o nes used in ISchneider fc Slusd (|2013h 
or in ISuvu et (l2014h . in order to assess the systematics 
due to the MST. However the principal problem remains: 
what is the possible range of “realistic” mass models? In 
other words, if a (possible composite) model k hts the data, 
how can one rule out that the corresponding model kx (with, 
say, 0.85 ^ A ^ 1.15) is not “realistic”? The various density 
profiles shown in Fig.3 indicate that real galaxies most likely 
exhibit a fairly large range of density prohles, which makes 
the classihcation of models prohles into “realistic” and less 
realistic ones fairly difhcult, if not all impossible. 

Thereby biased estimates of Hq are not fundamentally 
due to the power-law assumption of the lens model but are 
caused by the MSD itself, whose impact cannot be elimi¬ 
nated simply by studying the properties of material in the 
surrounding light cone, or through the adoption of more 
sophisticated lens models even when extra constraints are 
available from extended image conhgurations. In order to 
obtain a reliable estimate of Hq using gravitational lensing, 
one must hrst break the MSD. 

One possibility of doing so is to use the fact that MST 
also modihes image magnihcation /r —>■ p/A^. If the source 
luminosity, and thus the absolute magnihcation , is known A 
can b e hxed. For AGN as sources, as shown bv iBauer et al.l 
ll2012ll . the variability-luminosity relation can be used to 
estimate source luminosities but only in a statistical way 
with large scatter. Alt ernatively, if Type la supernovae are 
gravitation ally lensed llChornock et al.ll2013l : lOuimbv et al.l 
I2013l . l2014h . direct measurements of u would be possible due 
to their nature as approximate standard candles. In such 
cases, the MSD can also be broken and the multiply-lensed 
Type la supernovae become ideal laboratories to study the 
true l ens mass distribution as wel l as cosmological param¬ 
eters llK olatt fc Bartelmaiinl 1 19981 : lOeuri fc Kawang l2003l : 
iBolton fc Burled 2003h . 


A second way to break the MSD is to include in¬ 
dependent information about the lens mass distribution. 


e.g., from stellar dynam i cs (e.g.. iBarnabe fc Koo pmanj 
20071: Barnabe et al. 20091. I 2 OIII : iTreu fc Koopmaiia 20021 : 


Suvu et al.ll2010l . l2013l ~ 2014| j. The problem is that the cur¬ 


rent accuracy (typically 10%) in the velocity dispersion mea¬ 
surements may not yet be sufficient to break the MSD, not 
to mention the uncertainty due to the anisotropy of stellar 
orbits, which can lead to a system atic error on the slope at a 
level of ~ 5% (lAgnello et al. l2013l j . Further investigations in 
this regard were not included in this analysis but are greatly 
encouraged. 
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Figure Al. The slopes defined within different radial ranges are compared for the lens sample with = 0.6 and Zs = 1.5. On the 
left, the slopes of the local convergence k{ 6 ) measured between 9i = O.S^e s-nd 62 = 1-20e are plotted versus those measured between 
0 -^ = O.50E and 62 = I-S^e? on the right the relation of the cumulative distribution 9) is given. 


APPENDIX A: VARYING THE RANGE OF THE STRONG TENSING REGION 

Different angular ranges [^ 1 , O 2 ] will not result in exactly the same transformed power-law profiles and the associated As. It is 
therefore interesting to consider if our final statistical results strongly depend on the choice of the image range. In Fig. ixn we 
present the slopes measured between 9i — 0.50b and 02 = 1.50b and between 0i = 0.80b and 02 = 1.20e. The corresponding 
A distributions are shown in Fig. IA2I As can be seen, although measurements in these two cases are not exactly the same on a 
one-to-one basis, the statistical distributions of A are independent of the image range sets. We have, therefore, set 0i = O.50e 
and 02 = 1.50b to calculate the relevant quantities throughout the paper. 


APPENDIX B: DENSITY PROFILES 

In this appendix, we show the density profiles of our galaxy samples. Fig. IbT] presents the shape parameter distribution in 
the s — 5 plane (see Sect. 2 for definition). Fig. IB2I shows the surface density profiles of galaxies from each subgroup of the 
lens sample presented in Fig. 3, which are now normalized by galaxy effective radii 0eff. 


APPENDIX C: THE A DISTRIBUTIONS OF ALL LENS SAMPLES 

In this appendix, we present the distributions of A (and A) as a function of (tsib and of the “measured slope” s\ (and s\) for 
our lens samples with different sets of Zd and Zs combinations. 
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Figure A2. The distributions of A versus sx (top panels) and A versus sx (bottom panels) for the lens sample with = 0.6 and 
Zs = 1.5: the left-hand side panel shows the distributions calculated assuming 6 i = O.50e ^2 = foi* comparison, the results 

for which 9i = O.S^e s-iid ^2 = 1-20e are adopted are given on the right-hand side. On top of the scattered data (blue dots), the black 
squares (and red stars) with the error bars indicate the mean (and median) and its standard deviation within each bin. The dashed lines 
indicate where A = 1 and where A = 1; while the dotted lines indicate where sx = ^ and where sx = 1- 
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Figure Bl. The s — ^ distribution of the selected lensing galaxies at = 0.2, assuming = 1.5. Four colours represent four subgroups 
divided according to s and ^ of the local convergence distribution k{9): subgroup I (red) are defined as galaxies with s > 1.2, ^ > 1.02; 
subgroup II (blue) are those with s > 1.2, ^ < 0.98; subgroup III (orange) represents galaxies with s < 0.95, ^ < 0.98; and subgroup IV 
(green) are those with s < 0.95, ^ > 1.02. Contours indicate where the transformed slopes s\ = [0.3, 0.5, 0.7, 0.9, 1.1, 1.3, 1.5, 1.7]. 



0/0eff B/e,„ e/e,„ 


Figure B2. The surface density profiles of galaxies from each subgroup of the lens sample presented in Fig. 3 are shown as normalized by 
galaxy effective radii 0eff- The black curves show the total surface density distribution /c(0), scaled up by a factor of 5 for clarity; the red 
and blue curves represent profiles of the projected dark matter /^(^)/dm(^) projected baryonic mass k.(9){1 — /dm(^))? respectively; 
green dots indicate where the Einstein radii are with respect to the effective radii. 
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Figure Cl. Distributions of A (and A) versus (Jsie for six samples with and Zs indicated in the panels. On top of the scattered 
data (blue dots), the black squares (and red stars) with the error bars indicate the mean (and median) and its standard deviation within 
each bin. The dashed lines indicate where A = 1 and where A = 1; while the dotted lines indicate where = 1 and where sx = 1. 
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Figure C2. Distributions of A versus s\ and A versus s\ for all six samples with and Zs indicated in the panels. On top of the scattered 
data (blue dots) the black squares (and red stars) with the error bars indicate the mean (and median) and its standard deviation within 
each bin. The dashed lines indicate where A = 1 and where A = 1; while the dotted lines indicate where = 1 and where = 1. 
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